Effective medium theory of negative index composite metamaterials 
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In the homogenization of composite metamaterials the role played by the relative positions of 
the wires and resonators is not well understood, though essential. We present a general argument 
which shows that the homogenization of such metamaterials can be seen as the homogenization of 
ID single negative stacks. The sensitivity to the geometrical parameters is due to the fact that light 
sees the harmonic mean of the dielectric constant when propagating parallel to the layers. Since the 
dielectric constant is not positive definite the harmonic mean is not bounded and presents abrupt 
variations. We discuss applications of this remarkable phenomenon. 



Meta-materials with a possibly negative index of re- 
fraction have been under intense study since the results 
of Pendry 0, Q suggested that they might be fabricated 
by interspersing metallic wires (giving an effective nega- 
tive permittivity) with magnetic resonators such as split 
rings or dielectric fibers (giving an effective negative per- 
meability) . However the details of the exact way in which 
these two effects mix and perturb each other are not 
well understood. In this letter we present an effective 
medium theory that accounts quantitatively for this be- 
havior, provides an explicit geometrical recipe for obtain- 
ing true negative index media and opens the way to new 
applications. 

No general argument exists that by mixing magnetic 
resonators and wires one must obtain a negative index 
material at some frequency. Indeed, as the results of 
Pokrovsky and Efros have shown the negative ep- 
silon effect of metallic gratings is fragile with respect to 
the brutal mixing of the metallic wires with other ob- 
jects. Negative epsilon is essentially an interference phe- 
nomenon and anything that perturbs the interference of 
the waves scattered by neighboring wires is liable to de- 
stroy it. The wires need to "see" each other. The results 
of Pokrovsky and Efros are the first indication that the 
effective medium depends critically on the microscopic 
geometry Q|. They essentially tell us that resonators 
should not be placed in between the wires, or the ag- 
gregate may lose its negative permittivity. The regions 
occupied by the wires must have higher connectivity, and 
the simplest way of acheiving this is a stack of alternat- 
ing rows of wires and resonators. Each row of wires will 
then act as a homogeneous slab with negative permittiv- 
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row of resonators will act as a homogeneous slab with 
negative permeability (fx = 1 — j^ji , with a and u> mp 

determined by the geometry Q,@). We thereby arrive 
at the remarkable and quite general conclusion that the 
homogenization of complex wire/resonator composites is 
the homogenization of ID stacks of alternating single neg- 
ative homogeneous slabs. This conclusion is a direct con- 
sequence of the fact that the negative permittivity of wire 
gratings depends critically on their mutual interferences, 
in contrast to the negative permeability effect which is es- 
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Figure 1: The structure we study is a ID periodic medium 
with a unit cell composed of two layers, one with negative e 
(and positive /i), one with negative \i (and positive s). The 
H field is directed in the Z direction and the incident field is 
directed upwards. 



sentially a resonance phenomenon and does not depend 
critically on the coupling between neighboring resonators 

HQ 

However, while some interesting results on ID single 
negative media have been published 0, 0, 0] , the nat- 
ural question of the precise conditions under which such 
structures can be homogenized to obtain a negative index 
medium has not been satisfactorily answered. 

In this work we offer a rigorous self-contained homoge- 
nization analysis that gives a clear description of the ID 
alternating single negative medium described above and 
pictured schematically in Figure f . We shall use the term 
u P-k negative" to unambiguously refer to waves or homo- 
geneous materials that may support waves which have a 
dot product of the Poynting vector and wave vector that 
is negative. 

The discriminating factor, as in any homogenization 
problem, is the ratio between the wavelength and the 
period of the structure. In our case a complication is 
introduced by the extreme dispersion of resonator be- 
havior, but we avoid it by fixing the wavelength (and 
therefore the internal resonator geometry) and scaling 
the relative placement of wires and resonators. The dis- 
criminating geometrical parameter of the homogeneous 
model is therefore the ratio r = 4 where A is the free 
space wavelength and d is the period. This parameter 
is, a priori, large, which authorizes us to use the well- 
established two-scale expansion approach. In the follow- 
ing we treat only the H polarization case. The results 
hold in E polarization as well, provided we exchange fi 
and £ everywhere. 
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We note the macroscopic variables with capital letters 
X and Y and the microscopic variable with the small 
letter y as in Fig. 1. We first need to introduce the small 
parameter r\ — d/X which will be used in the expansion. 
Next we define the microscopic variable y = — which 
varies fast, ij being small, so it describes the system on 
the microscopic scale, the scale of an individual layer. 
Since we are interested in the homogeneous behavior of 
the system, at the end of the analysis this variable will 
be averaged over and will disappear. It can be regarded 
as an internal degree of freedom of the unit cell. 

For a ID structure with piecewise continuous e(Y) and 
Maxwell's equations can be written as 



1 dU 
dY \e(Y)dY 



k 2 sin 2 (6) 



U(Y) = 



e(Y) 

(1) 

where H(Y) = U(Y, y ) e ^sin(e)^ k = 2lr /x, 9 is the angle 
of incidence in the z = plane and where we will note 
fc 2 sin (0) — kx below. This second order differential 
equation can be written as a system of two first order 
equations 
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where we have introduced the new variable £ defined by 
equation and where special care must be taken with 
the derivative operator as it now takes the form 
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We develop the quantities of interest in multi-scale se- 
ries with parameter rj 



U(Y) = U (Y,y) 
£(F)=£o(^y)- 



-vUi(Y,y) 



(5) 



and insert them into equations El and El Note that all 
quantities on the right side of the equations must be pe- 
riodic in y with period d. By identifying same order 
terms in rj we obtain the following relations. From Eq.[2] 
we obtain 

~0 and e£ = _ + _ (6) 



while from Eq. El we obtain 



dy 



and 



dY 



dy 



M-f )r„ = u (7) 



We notice that the zeroth order terms of the U and £ 
expansions do not depend on y, the microscopic variable. 
However the first order terms do depend on it and are 
periodic of period d as noted above. This must be kept 



in mind as we proceed to average the equations (to ob- 
tain the homogeneous behavior of the fields) because the 
derivative of any continuous periodic function averages 
to zero. We note the average over y of a quantity / as 

(f) = ifof(y)ty- 

Taking the average ofEJandQthe y derivatives disap- 
pear as explained above and we get 

/ \ t dU ° 
{£Ho = ^Y 



dY 
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k x )U = 0. 
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which we put back together into one second order equa- 
tion 



d 2 u 

dY 2 



k 2 (s) M 



k 2 x \Uo = 0. 



The quantity in parentheses is the propagation con- 
stant in the Y direction, k\ and we thus obtain the dis- 
persion relation of the homogeneous medium: 
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where 



(/i) , e x = (e) , e Y = (e x ) 1 



(9) 



(10) 



The relations summarize the effective medium 
model and the rest of the letter is concerned with ex- 
ploring their consequences. The most important of them 
are: 1). The medium is anisotropic, 2). Its anisotropy 
can be extreme due to the unboundedness of the har- 
monic mean, 3). This behavior is strongly dependent on 
the microscopic geometry of the composite. We offer ex- 
plicit expressions illustrating the third point below while 
the first two are discussed after the presentation of the 
numerical results validating our model. We should note 
that we have not made any assumption about the losses 
so, in general, all quantities inE3 are complex. 

Eqs. El and El show that for H polarization, light sees 
the arithmetic mean of epsilon when propagating in the 
Y direction but the harmonic mean when propagating in 
the X direction. Let us denote the thicknesses of the two 
layers hi and h 2 (we assume hi + /12 = 1 for simplicity) 
and the corresponding permittivities £1 = ai + ibi and 
c2 = 02 + i&2- Rewriting Eq. E3we have: 



ex — hisi + /i2£2, ey 



eie-2 



hie 2 + h 2 ei 



(11) 



In this form it is easiest to see how the three main conse- 
quences come about. The first one is obtained by setting 
ex = £y- One quickly arrives at E\ = e 2 which contra- 
dicts our assumptions. The second one is seen by noting 
that there is nothing preventing the denominator of e Y 
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from approaching zero, ignoring losses. The third results 
from the prominent role played by the geometrical pa- 
rameters hi and hi. Their importance becomes evident 
as we explore the conditions under which the material is 
P ■ k negative. 

k 2 k 2 

Maxwell's equations tell us that P ■ k = H — — in 
a uniaxial anisotropic material in H polarization. Com- 
paring to Eq. El we obtain P ■ k = & 2 /i e ff- The suffi- 
cient condition to have a P ■ k negative material in H 
polarization is therefore that the effective permeability 
be negative. If in addition we require kx and ky to be 
real (ignoring absorption) for all angles of incidence we 
must also impose ex < and ey < 0. Writing it out 
explicitly using Eq. El we obtain the rather restrictive: 



a) 



hi 
hi 



< 



< 



< when e l < and < 

< £ when e 2 < and < 



(12) 



The crucial role played here by the purely geometrical 
quantity hi / hi combines with the unboundedness of ey 
noted above to explain the critical dependence of the be- 
havior of the composite on the relative spatial placement 
of its building blocks 0,0]. Moreover, these relations 
give, for the first time to our knowledge, a precise geomet- 
rical recipe for the conditions under which a composite 
will exhibit true negative index properties. We believe 
that, along with Eqs. El from which they are derived, 
they will become useful in the theoretical and especially 
the experimental study of negative index composites. 

We now study the effect of losses on the effective 
medium. We assume the media are passive (b > 0), 
and that the imaginary parts are comparatively small 
(b <C \a\) so that we can ignore quadratic terms in b such 
as 6162, etc. Rearranging El we get 

ex = (a) +i {b) 

1 -( ba ~ 2 ) 
e Y = (a l ) +i\^r- 

Let us take a closer look at ey. The real part is the 
harmonic mean of the real parts of the permittivities of 
the layers. These permittivities have different signs, and 
therefore this quantity is not bounded a priori. Let us 
call it A. Rewriting ey we get 



ey 



A 



(13) 



The unboundedness of the harmonic mean has a much 
stronger impact on the imaginary part than on the real 
part. Eq. El is only valid when max(6i,62) <C A -1 . 
Other cases are of lesser interest as losses dominate. 

We now validate our model as formulated in Eq. El 
with numerical simulations. The transmission of the ho- 
mogeneous anisotropic slab was computed using a ID 
transfer (T) matrix method, while for the alternating 
single negative stack we chose a scattering (S) matrix 
method for its numerical stability. In what follows all an- 
gles and phase differences are measured in degrees while 
the scale parameter r is dimensionless. 
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Figure 2: Offset between the power transmission coefficients 
of single- negative stack (||7a|| 2 ) and the equivalent homoge- 
neous medium (||T£r|| 2 ) as a function of the scale parameter 
r = X/d for normal incidence (a) and of angle of incidence for 
r = 35 (b). They are within 0.05 of one another for r > 30 
for all angles of incidence, c) Phase mismatch between the 
transmission coefficients of single negative stack and effective 
medium as a function of incidence angle for r — 35. They are 
within 7° of each other for all angles of incidence. 



The structure under study is composed of 200 alter- 
nating layers of thicknesses hi = 0.4, hi = 0.6, with 
permittivities £i = 2, si — —2.8 and permeabilities 
Hi = —3.295, fi2 — 0.53. The equivalent parameters 
resulting from Eqs. El are £ x = —0.4, ey = —70 and 
H e ff = —1. The total thickness of the structure is 100 
periods. We now present numerical simulations com- 
paring the transmission coefficient of the ID alternating 
layer structure, Ta, to that of the equivalent homoge- 
neous structure, Tjj, as a function of r and of angle of 
incidence. 

Figure 2. a) shows the offset between the power trans- 
mission coefficients, ||Tif|| 2 — ||Ta|| 2 , as a function of the 
scale parameter r for normal incidence. The same off- 
set as a function of angle of incidence appears in Fig. 
2.b) for r — 35. Finally, in Fig. 2.c) we plot the phase 
mismatch between the two transmission coefficients as a 
function of angle of incidence for r — 35. These results 
demonstrate clearly that the meta-material behaves in 
all respects as a homogeneous medium with the effective 
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parameters given by Eqs.E3f° r periods more than 20 to 
30 times smaller than the wavelength. 

The behavior at shorter wavelengths can be described 
by pushing the developments of Eq. Qj]to higher orders. 
The simple effective medium description is then no longer 
quantitatively valid but, interestingly, negative index- like 
behaviors remain More precise results in this direc- 
tion will be presented elsewhere. 

We have this far detailed only the last of the three main 
consequences of Eq. It is expressed quantitatively 
in Eqs El and El We now discuss the first two: that 
negative index resonator/ wire composites are anisotropic 
and that their anisotropy can become extreme due to the 
divergence of the harmonic mean. 

The significance of the anisotropy of the effective 
medium resides in the difficulty of interpreting experi- 
mental data from a single angle of incidence. The method 
of extracting effective parameters that is most commonly 
used in the literature [lj, LLJ, LLj, LLg] may fail when faced 
with an anisotropic material, particularly when it is in 
addition indefinite. This aspect points to the need for a 
new extraction procedure, one that can unambiguously 
place the sample under study within one of the cases 
listed by Smith and Schurig ^jj a t a given frequency. 
Such a procedure would necessarily involve transmission 
and phase measurements over multiple angles of inci- 
dence. 

The implications of the second point stretch beyond 
the anisotropy of resonator/ wire composites. Among 
other things, it tells us that it is possible to realize mate- 
rials with very large positive or negative permittivities or 
permeabilities using low valued single negative stacks. It 
also enables a new understanding of enhanced transmis- 
sion through very thin slits in thick metal screens by 



treating them as ID air/metal stacks, providing a gen- 
eral theoretical framework for the results of Shen et. al 
[lflj . Another promising application is related to plas- 
mon lasing (2flEll |. Our preliminary results show that by 
stacking thin layers of active material with thick metal- 
lic layers it may be possible to use the divergence of the 
harmonic mean to dramatically enhance the gain of the 
active material. More precise results and examples will 
be published elsewhere. 

In conclusion, we have argued that, because of the 
physical nature of the negative permittivity phenomenon 
in metallic gratings, composites employing them may be 
simply modeled as ID single negative stacks. Rows of 
wires make up effective negative permittivity slabs while 
rows of resonators (metallic or dielectric) provide the 
negative permeability. We have shown using a rigorous 
self-contained homogenization theory that light sees the 
harmonic mean of the dielectric constant when propa- 
gating parallel to the layers. Since the dielectric con- 
stant is not positive definite the harmonic mean is not 
bounded. This remarkable phenomenon, never before 
pointed out, explains quantitatively the sensitivity of the 
composites' behavior to their microscopic geometry and 
provides clear criteria for constructing a true negative in- 
dex metamaterial. Furthermore, it tells us that it may 
now be possible to fabricate media with arbitrarily large 
positive or negative effective parameters using conceptu- 
ally simple and technologically straightforward ID sin- 
gle negative stacks. More in depth investigations of the 
possible applications of this phenomenon are promising 
directions for further research. 

This work was supported by the French RNRT project 
CRISTEL and by the EU-IST project PHAT 510162. 
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